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Abstract. We further investigate the long-term evolution 
of a trajectory of a stellar-mass orbiter which is gravita- 
tionally bound to a massive central body acting by Newto- 
nian force. The orbiter undergoes repetitionary collisions 
with an accretion disc. We consider eccentric orbits inter- 
secting the disc once or twice per each revolution and we 
solve the equations for osculating elements. We find the 
terminal radii of the orbits and time needed to bring the 
orbit into the disc plane as function of initial parameters, 
and we show that previous simplified estimates (derived 
for the case of low eccentricity) remain valid within factor 
of two. The discussion presented in this paper off'ers a toy 
model of the orbital evolution of a satellite which passes 
through the rarefied gaseous environment of a stationary 
disc. We demonstrate that the drag on the satellite should 
be taken into account in calculations of the stellar distribu- 
tion near the super-massive black hole in galactic centers. 
On the other hand, it does not impose serious restrictions 
on the results of gravitational wave experiments, except 
in the case of discs with rather high surface density. 
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1. Introduction 

Rapid stellar motion in galactic centers reflects properties 
of a presumed supermassive compact body in the nucleus 
(Binney & Tremaine 1987; Merrit & Quinlan 1998a, b). 
One can define the sphere of influence, rh = G'M/(cr)^, 
which determines the region of space where gravity of the 
centre dominates over more remote parts of the galaxy 
(M ^ 1O^M0 is the central mass, cr « 10^ km/s is veloc- 
ity dispersion). Inside rh, a star can be considered as a 
satellite orbiting around the central body. We introduce 
gravitational radius rg = 2GM/c^, which characterizes 
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compactness of that body and provides a natural length- 
scale for the problem discussed in this paper. When ex- 
pressed in terms of rh and rg, we will explore in this paper 
the region rg ^ r ^ rh. 

Several mechanisms act as perturbations on almost 
free motion of orbiters in the central gravitational field 
inside rh: gravity of the background galaxy, dynamical 
friction in the field of other orbiters, energy losses due 
to gravitational radiation, tidal interactions, and friction 
caused by the rarefied gaseous environment (see Rees 1998 
for a recent review) . Here we will concentrate on the latter 
effect, namely, we will consider the long-term evolution of 
orbital parameters (eccentricity e, semimajor axis a, incli- 
nation i, argument of pericentre w) of the orbiter which 
crosses the plane of an accretion disc (Fig. |l|) . We as- 
sume that the disc is geometrically thin (collisions will be 
treated within the impulsive approximation) and neglect 
all other dissipative effects, ablation of the orbiter, correc- 
tions due to general relativity, as well as self-gravity of the 
disc. Some of these effects were examined quite recently 
by Vokrouhlicky & Karas (1998b), Ivanov et al. (1998), 
and Zurek et al. (1994). In this context, gradual evolu- 
tion of the disc surface density was explored by Ivanov et 
al. (1999). Further citations to older works can be found 
in Syer et al. (1991), Vokrouhhcky & Karas (1993), and 
Artymowicz (1994). 

The impulsive approximation ignores effects of reso- 
nances, and it assumes nonzero inclination of the orbit 
with respect to the disc plane. Our calculations thus do 
not apply to the case when the orbiter remains embedded 
inside the disc, the situation which has been widely dis- 
cussed because of its profound implications for describing 
the formation of stars and planets (Lin et al. 1999; Bryden 
et al. 1999). Also, we do not consider possible disruption 
of the orbiter by the central body; this process has been 
discussed by many authors, e.g. Kim et al. (1999) very re- 
cently. Under these assumptions, the orbital evolution can 
be solved in two steps: First, parametric relations for os- 
culating elements are obtained which determine the shape 
of the orbit (e(i), a(i), etc). Next, the temporal depen- 
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Fig. 1. Geometry of tlie model. The satellite (mass M,^Mq) 
follows a slowly evolving elliptical orbit around more massive 
central body {M~10^Mq). Argument of pericentre lo and incli- 
nation angle i specify orientation of the ellipse. In the present 
model, action of the disc matter upon the satellite depends 
only on the disc properties at the places of intersections. For 
Lu — -K 12 (as in this sketch) both intersections with the disc 
occur at the same distance, while the case is just opposite for 
oj = (the places of intersections then correspond to the peri- 
centre and the apocentre, respectively). 



dences are examined. Although the present discussion is 
very much simplified, it captures basic properties of the 
mechanism driving the orbital evolution, and it offers a 
toy model which can be used as a test bed for more com- 
plicated and astrophysically realistic studies. 

In this paper we extend discussion of Vokrouhlicky & 
Karas (1998a; cited as VK hereafter) who studied the or- 
bital evolution under the assumption that orbiter's tra- 
jectory intersects the disc only once per revolution. Such 
an assumption applies to the case when the disc has a 
finite radial extent and the orbiter follows an eccentric 
path (the initial stage of the capture from an energet- 
ically unbound state). Regarding the description of the 
disc matter and its mutual interaction with the orbiter, 
VK adopted several assumptions which are quite plausi- 
ble for this type of the toy model: a geometrically thin disc 
rotating at Keplerian velocity in purely azimuthal direc- 
tion. It was further assumed that the disc itself does not 
change with time; in particular, dissipative waves induced 
by the orbiter and the possibility of opening of a gap was 
not not taken into account. These simplifying assumptions 
are maintained also in the present paper. 

Gradual circularization and precession lead eventually 
to the orbit with two intersections per revolution; at this 
moment the analytical solution from VK ceases to be 
valid. Terminal stages of the evolution were treated by 
Ranch (1995) under approximation of zero and moderate 
eccentricity, but orbits with large eccentricity could not 



be fully examined in terms of an analytic solution, be- 
cause the two intersections can occur at different radial 
distances. One thus needs to specify surface density of the 
disc as a function of radius. We write down corresponding 
evolutionary equations and solve them numerically under 
the assumption of the power-law surface density profile. 
We examine the temporal evolution of orbital parameters, 
showing that a very special orientation of the trajectory, 
(jj = 7r/2 (when both intersections are at the same r), 
represents a stable configuration which can be treated an- 
alytically. In other words, our present discussion extends 
previous discussion of the initial period of the orbital evo- 
lution, and it concerns also the intermediate stage, when 
two intersections per revolution have already developed 
while inclination is still nonzero. 

2. Equations for orbital parameters 

The basic formulation of the problem follows VK, in par- 
ticular, the osculating elements of the orbiter are updated 
at every intersection with an axially symmetric disc, which 
is described by surface density Sd oc r* (s = const) and 
by orbital velocity v^ = v^ e^. It appears natural to con- 
sider Keplerian rotation of the disc, Vcf, = v^^ — y^GM/r, 
which fits with the assumption about its small geometri- 
cal thickness, and enables us to introduce vertically inte- 
grated quantities in the standard manner (we examined 
also different rotation laws; Karas & Subr 1999). The cen- 
tral body attracts the orbiter by Newtonian force. Pas- 
sages through the disc are treated as inelastic collisions 
which, by assumption, do not affect the properties of the 
disc and of the orbiter itself. Stated in a different way, disc 
material located along the orbiter's path transfers part of 
its momentum to the orbiter, but then the disc restores 
its original stationary state before the next collision oc- 
curs. Collisions thus lead to tiny changes of the orbiter's 
velocity, v ^ v' . The momentum conservation law gives 
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Vr Br + V^ 6^ + {v^ + AVj^) B^ 
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where A{r) = Sd^roiSj^w^ . Here, v^d is the relative 
speed of the orbiter and the disc matter, while S» is 
surface density ascribed to the orbiter and defined by 
E* = M^/ (ttR^) (quantities denoted by asterisk refer to 
the orbiter). At every collision, a tiny kick changes velocity 
according to Eq. (n|) and leads to corresponding changes 
of orbiter's specific (per unit mass) binding energy E, spe- 
cific angular momentum L, its projection Lz onto rotation 
axis, and radial velocity component Vr'. 
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here we introduced variables x — cosi, y — 1 — e^ , z ~ 
1 + e cos Lj, and we expanded right-hand sides of Eqs. (2)- 
(5) to the first order in A (since we deal with a small 
perturbation of the orbit). The collisions occur, in general, 
twice per each revolution of the orbiter, at the ascending 
and descending nodes, which will be denoted by subscripts 
"1" and "2" , respectively. Hydrodynamical drag acts when 
the orbiter passes across the nodes. By adding the two 
contributions one obtains 
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Eqs. (6)-(9) determine the orbital elements in a paramet- 
ric form, and we can adopt a; as a convenient parameter. 
Notice that the time dependence has not yet been involved 
explicitly. Temporal evolution can be examined by intro- 
ducing the orbital period At = 2t: a?" I "^ J^TTM . 

In further equations we take advantage of the fact that 
the assumed density profile, Ed = K (r/rg)'' E,, enables 
us relate A's at two different radii: 
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Initially, the orbit is energetically bound or parabolic 
(E > 0) with a = ao, i — io, and uj ~ ujq. Subsequent evo- 
lution decreases inclination and eccentricity, and it causes 
slow precession in uj. Naturally, there is a lower limit on 
inclination, which determines a moment when integration 
of our evolutionary equations must be terminated for at 
least two reasons. First, the nature of the problem be- 
comes different once the orbiter is fully embedded in the 
disc (Lin et al. 1999; Goldreich & Tremaine 1980). Second, 
the assumption (|l|) about impulsive collisions requires the 
relative speed to be supersonic. 

Regarding the latter condition, we can write, for 
a satellite on a quasi-circular Keplerian orbit, v^ ~ 
2.1 X 10i°(r/rg)"^/^cms-^ Therefore, i;,oi ~ 4.2 x 
10^°sin(i/2) (r/rg)~ ' cms^^, which is to be compared 
with the speed of sound Vs- For example, in the standard 



Shakura-Sunyaev accretion disc (Frank et al. 1992) with 
the central mass M = 10^ Mq, the accretion rate M = 
O.OlMsdd (AfEdd is the accretion rate corresponding to 
Eddington luminosity), and viscosity parameter a = 0.1, 
one obtains Vg ~ 1.27 x 10^ (''/^g) cms~^, and the con- 
dition Vici^Vg then reads sin(i/2) ::^ 3.1 x 10^"* (r/rg) ' . 
Inclined orbits can be further distinguished into two 
sub-cases: Those with (i) a single intersection with the 
disc plane per revolution, and (ii) two intersections per 
revolution. 

2.1. The initial phase: single intersection 

Before describing the results of numerical integration we 
recall that the parametric analytical solution was given in 
VK under the assumption that collisions occur only once 
per revolution, at the inner orbital node (while the more 
remote node was assumed to lie beyond the outer edge of 
the disc). This situation corresponds to the initial stage 
with large eccentricity and appropriate orientation of the 
ellipse in terms of uj. In this case, one can verify that Eqs. 
(6)-(9) with ^2 = are equivalent to Eqs. (20)-(23) of 
VK. The corresponding solution reads 
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where ro, Si, -82 are constants of integration. We recall 
that, in Eqs. (12)-(14), x stands as parameter. Time does 
not enter in mutual relations among a(x), y{x), and z{x). 
In accordance with the previous definition, x has a phys- 
ical meaning of cosine of inclination, and one can find its 
temporal dependence by solving implicit relation 
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dx. 
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2.2. The subsequent phase: two intersections 



The orbit eventually develops two intersections per rev- 
olution, and one has to resort to a numerical treatment 
of the problem. However, we checked (see below) for the 
wide range of starting parameters that the results do not 
differ substantially from the case with special orientation, 
<^o = 7''/2 (i.e. 2 = 1), which can still be treated analyti- 
cally. In addition, one can verify that this particular ori- 
entation of the orbit is linearly stable, or, in other words, 
the orbital evolution restores the initial argument of the 
nodes if slightly perturbed from to = tt/2. This conclu- 
sion holds for Ed decreasing with radius (s < 0), but it 
requires a further remark: It was shown (Vokrouhlicky & 
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Fig. 2. Semimajor axis o (ordinate, units of ao) versus incli- 
nation i (abscissa). Inclination decreases from its initial value 
io = 7r/4 down to zero, while a decreases from a = ao to the 
final a = flf which depends on position of the outer edge of 
the disc. First, the evolution follows the thick solid line S (sin- 
gle intersection per each orbital revolution), then one of thin 
solid lines (two intersections per revolution). Dashed curves in- 
dicate lines of constant time (arbitrary units). Initial tJo = 0, 
eo = 0.5. See the text for details. 



Karas 1998b) that the claim about the stability is not true 
when the disc gravity is taken into account. In that case 
the stability could be maintained only if the orbiter-disc 
interaction is substantially stronger than it follows from 
Eq. (y), perhaps due to turbulence or different nature of 
the collisions. In fact, one should take the disc gravity into 
account also for self-consistency reasons (although we ig- 
nore it here) , because both gravitational effects of the disc 
and its dissipative influence on the passing orbiter arise 
from the same distribution of matter. 
The solution of the w = 7r/2 case is 



Ci(l + x) 
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(16) 

(17) 
(18) 
(19) 



where constants Ci and C2 are to be determined from 
initial values oi x = xq and y = yo- 

For e = one obtains C2 = 0, and Eqs. (16)-(18) thus 
reproduce Ranch's (1995) formula (7) for the final value of 
semimajor axis of the circular orbit as function of its initial 
value and initial inclination. In the present notation. 



Fig. 3. Eccentricity e versus inclination i for the same initial 
parameters as in Fig. m. Here, jq ~ ''r/4, loq ~ 0. Again, the 
thick solid line shows the initial phase of single intersection, 
while thin lines correspond to the subsequent period with two 
intersections. 
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Fig. 4. Semimajor axis a versus inclination i (solid lines) for 
orbits with different io (initial inclination). Contours of con- 
stant time are also shown (dashed; time is given with each 
curve in arbitrary units). Other parameters as in Fig. 0. Tra- 
jectories in this figure start with a = ao, i = io, and circularize 
at a = iflo. Then they proceed, already with e = 0, down to 
terminal radius af at i = 0. 



Denoting the initial semimajor axis by ao — a{io), one 
obtains for the constant 



Ci = flQ (1 + cos io) = 4 ao cos'' (io/2) . 



(21) 



Temporal evolution (19) can be integrated after spec- 
ifying the surface density profile by parameter s: 



t = Tri^-^cf^""/, 



(7-4s)/2i 



(22) 



a = Ci (1 + a;) 



where I-y{x) —J (1 -I- a;) '^^/l — xdx. Integration can be 
(20) carried out in terms of elementary functions for suitable 
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values of s (in particular for s = —3/4) by applying recur- 
rent formula 2(7— 1)/^ = (7 — 3)17-1 — {1 + xY' 



'VT 



and Ji/2 = arcsinx. For other values of s, numerical inte- 
gration is straightforward. 

3. Results 

We solved Eqs. (6)-(9) with different radial profiles of the 
surface density and different rotation laws of the disc mat- 
ter. Numerical analysis provides us with estimates of ac- 
curacy of the previous results, which were restricted either 
by assuming low eccentricity (the final evolutionary stage 
of the orbit), or, in turn, by ignoring collisions at the or- 
bital node which is farther away from the centre (the ini- 
tial stage). In other words, we relaxed some assumptions 
about eccentricity and orientation of the orbiter's trajec- 
tory 

Fig. shows semimajor axis as a function of inclina- 
tion, a{i). The orbital evolution proceeds from the initial 
state, ao = a{io), to the final af at i = (when orbiter's 
orbit plunges into the disc plane). Thick solid line (de- 
noted by "S" ) corresponds to the evolutionary stage with 
the single intersection, which is governed by Eqs. (12)- 
(14). (We checked that this analytical solution, valid dur- 
ing the phase of single intersections, is reproduced also by 
the numerical solution which will be used further.) The 
orbital evolution follows the line S to the point when two 
intersections occur. The actual moment when this hap- 
pens is determined by initial orientation of the ellipse and 
by the radial extent of the disc. Subsequent evolution is 
disjointed from the thick line, and it goes along one of 
thin solid lines, which were constructed according to Eqs. 
(16)-(18). Only this latter part of the solution depends on 
the value of s (here we set s = —3/4, corresponding to the 
outer region in the standard thin disc model). Dashed lines 
correspond to t = const in Eq. (19). They are plotted with 
constant time steps of At = 2 (arbitrary units) as indi- 
cated in the graph. For clarity of the graph, we do not show 
the lines of constant time near i = because evolution of 
the system slows down visibly when inclination is small. 
Notice that the final radius of inclined orbits ranges be- 
tween 0.48 < Uf/ao < 0.57, lower values corresponding to 
the case when the branches with two intersections develop 
late in the course of evolution. One can conclude that even 
if the moment when two intersections per revolution occur 
was unknown (e.g. because location of the outer disc edge 
is left unspecified) , still the final radius could be estimated 
with the relative error less than about 0.15. In this plot 
we choose loq = 0, because such initial orientation of the 
ellipse captures a long period of evolution when the ellipse 
has a single intersection per revolution. We recall that for 
a; = the intersection is at the pericentre, while the other 
passage across the disc plane occurs at the apocentre; cf. 
Fig. 0. 

Fig. |3| is constructed in analogous way as Fig. g, show- 
ing eccentricity e as function of i. All the curves of e{i) 



eventually meet together at the origin, which corresponds 
to circular orbits of zero inclination (but different af's). We 
constructed similar graphs as those in Figs. ||-^but for dif- 
ferent starting parameters, including the counter-rotating 
orbits; the results of such analysis are qualitatively the 
same and can be quantified by Of . Complementary to Fig. 
H, Fig. Q shows a{i) for different zq. Again, for the sake of 
definiteness, the initial ellipse has cq = 0.5, loq = 0, and 
the initial pericentre is identified with the outer edge of 
the disc (i.e. radius above which the disc density is con- 
sidered as negligible) . Such orientation of the orbit could 
be called as fiducial in our examples concerning with the 
capture of a satellite from a remote eccentric trajectory, 
initially extending well above rh. All the orbits in this ex- 
ample have a single intersection, which is located exactly 
at the outer disc edge until a decreases down to half of its 
initial value, a = iag (denoted by the horizontal dotted 
line). At this point the second intersection occurs for the 
first time. 

Terminal radius Of of a general orbit, with either one 
or two intersections per revolution, follows from the nu- 
merical solution, as it was described in previous para- 
graphs. Eqs. (16)-(17) represent a good approximation 
to the numerical solution of Eqs. (12)-(14), which be- 
comes exact for cjo = 7r/2. With this orientation of the 
trajectory we recover the above-quoted formula, origi- 
nally derived for moderate eccentricity (Ranch 1995): 
af~-jao (1 — Cg) (1 -t- cosio) . A simple formula has obvi- 
ous practical advantages, and one thus wants to quantify 
its accuracy. We define the relative error as 

, , . ^ ai{aQ,eo,iQ,ujo) - a{{ao,eo,io,TT/2) 
Aaf(eo,«o,^o) = (23) 

and show the surface plot of Acf in Fig. g for iq = 5°, 
s — —3/4, and different starting loq, cq. It is evident from 
the plot that the difference is substantially reduced when 
loq — > 7r/2 and/or 69 — » 0. Although this plot was con- 
structed for the particular value of ip, it turns out that 
the difference is also reduced for higher initial inclination. 
Our discussion is completed by calculating temporal 
dependences. Time evolution of semimajor axis is shown 
in Figs. ^0 for the disc with Sd oc r^^/^. Here we choose 
luq — 7r/2, so that the orbit has two branches from the very 
beginning. Temporal dependences can be computed by di- 
rect evaluation of Eq. (19), or Eq. (|2|) when it applies, but 
purely numerical integration of original Eqs. (6)-(9) gives 
the same results. Characteristic time-scale r (in terms of 
orbital periods) is determined by ratio of surface densities, 
T — S]*/I]d. When scaled to the form appropriate for the 
gas-dominated region of the Shakura-Sunyaev disc around 
a supermassive compact core (Frank et al. 1992), 
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Fig. 5. Relative difference Aaf(eo,(^o) of final radius, as de- 
fined in Eq. (psl). This figure illustrates the fact that Eqs. 
(16)-(19) estimate the correct value of af with accuracy better 
than 15% in the whole range of eccentricity. 



Fig. 7. Time evolution of semimajor axis as in Fig. |6| but 
for io — it/4: and different initial eccentricities (solid lines; eo 
is given with each curve). Isolines of constant inclination are 
shown with constant step 0.05 (dashed). 




Fig. 6. Time evolution of semimajor axis a{t)/ao (solid lines). 
Initial inclination is given with each curve. Isolines of constant 
eccentricity (dashed) indicate how the pace of orbital evolution 
eventually slows down in the case of orbits which are almost 
circular and have small inclination. Again, time is in arbitrary 
units which can be scaled with the disc surface density (see the 
text); eo = 0.5, loo = 7r/2. 
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From here, the fiducial value of t corresponding to a solar- 
size satellite (E^wA/Q/rg) passing across the inner part 
of disc (where density has maximum) is of the order 10^ 
revolutions (inverse of the numerical factor in front of rhs. 
of Eq. (24)). Scaling with the model parameters is also ev- 
ident: e.g. for an orbiter at r « lO^^g, and the disc with 
viscosity parameter a « 10~^, accretion rate M « IMq 
per year, E* = Eq, one obtains r « 10^ revolutions, corre- 
sponding to w 10* yrs. Although the precise time-scale can 



be computed in a straightforward manner from Eq. (19) 
above, assuming the disc model, one should be cautious 
about interpretation of time intervals in physical units. 
In a realistic model, major uncertainty affecting the con- 
stant of proportionality K is connected with poorly known 
effects of turbulence (Zurek et al. 1996), which tend to 
enhance the drag on the orbiter. We recall that K stands 
only in front of the time integral, and it determines tempo- 
ral dependences in absolute units, but it does not appear 
in mutual relations among orbital parameters. Knowledge 
of K is thus unnecessary for determination of the shape 
of the orbit and its evolution in arbitrary units of time, 
provided that collisions can be treated as small perturba- 
tion. 

In the above-given discussion we restricted ourselves 
to Keplerian discs. This turns out to be rather important 
assumption because non-Keplerian discs can transfer en- 
ergy to the orbiter. In such a case, a(t) passes through a 
minimum at some point and starts increasing again. For 
example, such behaviour occurs in a disc with constant 
angular momentum density / (Karas & Subr 1999), if one 
can still retain other assumptions of our model, namely the 
impulsive approximation for collisions (/ = const appears 
in discussion of discs which acquire substantial geometri- 
cal thickness, so that vertically integrated quantities may 
be too inaccurate). 

4. Conclusions 

There are two areas where the models, like the one which 
we discussed here, are relevant. First, the calculations of 
the stellar distribution in galactic cores (Bahcall & Wolf 
1976; Quinlan et al. 1995) predict the presence of density 
cusps and corresponding integrated-light profiles around 
the massive black hole. Stellar orbits and, consequently. 
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the overall form of the cusps will be modified by the pres- 
ence of gaseous environment. We demonstrated that evo- 
lutionary time-scales are typically 1-2 orders of magnitude 
less than the stellar life time, and the influence of the disc 
should be thus accounted. This is true especially for ac- 
tive galactic nuclei with enhanced accretion rates, while in 
the center of Milky Way the gaseous environment is too 
rarefied to exert any substantial hydrodynamical drag on 
stars. The smallest resolved scales (« 0.1") do not allow 
direct observational comparisons at present, but this may 
be possible when better resolution is available (although 
the direct tracking of individual stars does not seem pos- 
sible in near future; Kormendy & Richstone 1995). 

Another possible application concerns the predicted 
gravitational waveform patterns from inspiralling satel- 
lites. It was suggested (Chakrabarti 1996) that the wave 
emission may be significantly altered if the satellite moves 
inside a super-Eddington (M « lO'^AfEdd) flow and passes 
through shocks there. This could be an interesting possi- 
bility but it represents rather exceptional situation. When 
most of the motion takes place outside the disc, as in 
our case, the time-scales of the orbital evolution due to 
the presence of the disc are always much longer than 
the orbital period, even for quite high accretion rates 
(M Ri AfEdd), and the production of gravitational waves 
is virtually unaffected. One should however notice that 
gaseous environment helps to capture the satellite and 
bring it down to the low orbit with small inclination, the 
situation which is needed for the efficient wave emission. 

We solved Eqs. (6)-(9), which determine the orbital 
evolution of a body on an eccentric trajectory evolving 
adiabatically due to dissipative collisions with the disc. 
In the general case it is quite obvious from the form of 
these equations that one has to resort to the numerical 
approach, but we checked with discs obeying reasonable 
power-law surface density profiles (—1 < s < —0.5) and 
Keplerian rotation law that evolution of the true trajec- 
tory is well approximated by the corresponding orbit with 
oJo — 7r/2. This case can be integrated analytically [Eqs. 
(16)~(19)]. Approximate terminal radii of circularized or- 
bits agree with those obtained by numerical integration 
within precision better than 15%. This conclusion can be 
put in other words, namely, effects of the disc self-gravity, 
of different rotation laws, and of mutual collisions among 
orbiters seem to infiuence the evolution more substantially 
than do exact properties of the disc material. 

Compared with previous papers, our present discussion 
has been more general on the level of semi-analytical treat- 
ment of the orbits with arbitrary orientation, for which 
either one or two intersections develop in the course of 
evolution. In particular we examined how the power-law 
index of the disc surface-density profile enters in equations 
for the orbital shape and time evolution of the parameters. 
It should be quite obvious that we picked up one aspect 
of a complex problem of the long-term orbital evolution, 
which will be build into a mosaic to form the complete 



picture in astrophysically more realistic models. Various 
compelling processes were already considered separately. 
Now we need to take them jointly into account, especially 
the effects of resonances between the satellite and the disc 
matter, corrections due to general relativity (of the gravi- 
tational field of the central object and energy losses of the 
orbiter), tidal effects, and gravitational interaction with 
other satellites. This will certainly require mostly numer- 
ical approaches for which analytical estimates of individ- 
ual processes and toy models, like the present one, provide 
useful guidelines. 
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